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ABSTRACT: We apply here a recently developed approach to compute 
the short distance corrections to scaling for the correlators of all primary 
operators of the critical two dimensional Ising model in a magnetic field. 
The essence of the method is the fact that if one deals with O.P.E. Wilson 
coefficients instead of correlators, all order I.R. safe formulas can be obtained 
for the perturbative expansion with respect to magnetic field. This approach 
yields in a natural way the expected fractional powers of the magnetic field, 
that are clearly absent in the naive perturbative expression for correlators. 
The technique of the Mellin transform have been used to compute the I.R. 
behavior of the regularized integrals. As a corollary of our results, by com- 
paring the existing numerical data for the lattice model we give an estimate 
of the Vacuum Expectation Value of the energy operator, left unfixed by 
usual nonperturbative approaches (Thermodynamic Bethe Ansatz). 
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1 Introduction 



It is well known that (D = 2) conformal field theories, ffl, can be used to 
describe statistical systems at the critical point, @. The second step from 
this point of view is to obtain informations on the behavior of the system near 
criticality, that is to study a conformal field theory perturbed by a relevant 
operator (i.e. with scale dimension < x < 2). 

On this line a wide class of perturbations have been found to give rise to 
integrable models, ||, i.e. models in which an infinite number of conserved 
charges exists and constrains the S matrix to be factorized (and possibly elas- 
tic). Starting from the exact knowledge of the S matrix one can investigate 
the off shell behavior of the theory by use of the form factor method |4j] that 
essentially gives rise to a long distance expansion for correlators in terms of 
er MlT , Mi being the masses of the complete theory and r the argument of 
correlators. 

Another approach to the problem [^, ^[ is obtained building up in some 
way a perturbation theory around the (massless) conformal field theory, from 
which one could get informations on the short distance behavior of the com- 
plete theory (the small adimensional parameter being now the coupling con- 
stant times an adequate positive power of r). This line of research is im- 
portant not only because it is complementary to the previous one, but also 
because of its generality (integrability is not essential, thus all perturbations 
can be treated if the starting point conformal theory is known) . However it is 
well known that when the coupling is relevant, naive perturbative expressions 
for correlators are plagued by I.R. divergences, that constituted an unsur- 
mountable obstacle to the construction of a general theory of perturbations 
in these cases. Nevertheless in f7| some recent results results f|, ||, [l(| (see 
also [|11], 0, [13|] for preexisting ideas) have been developed to give an all order 
I.R. safe general approach to describe the short distance behavior of confor- 
mal theories perturbed by relevant operators. The main idea of the method 
is the fact that Wilson coefficients of Operator Product Expansions, being 
short distance objects, can be taken to have a regular, I.R. safe, perturbative 
expansion with respect to the coupling. 

In this paper we apply the above mentioned OPE technique to reconstruct 
the short distance behavior of the (continuous) critical two dimensional Ising 
model perturbed by a magnetic field H: clearly our results can be applied to 
describe the scaling limit of the corresponding model on the lattice. 

We will derive the first nontrivial terms of the expansion in the scaling 
variable t oc HR 15 ^ 8 (R, H being respectively the lattice spacing and lattice 
magnetic field) for all the correlators of the primary operators of the critical 
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conformal theory. In particular we obtain naturally the predicted [fT4 , 
^8/15 correc tions to the spin spin correlator. 

The OPE approach leaves in general unfixed some constants that parame- 
terize the Vacuum Expectation Values of the used low dimensional operators. 
One of these constants has been fixed by use of known non perturbative re- 
sults (Thermodynamic Bethe Ansatz), while the other has been estimated by 
comparison with the available data from numerical Monte Carlo simulations 
on the lattice. 

It is worth remarking here that the perturbative expansions in H/\T — 
T c \ 15 ^ 8 , (obtained starting from the massive noncritical Ising at zero field) 



that are presented in are clearly inapplicable when T = T c : this is 
another manifestation of the I.R. divergences we spoke above. We notice 
also that, being the critical Ising Model in magnetic field an integrable model 
with known S matrix, its long distance behavior can be computed by use of 
form factor techniques ITB, 17 ]. 



The plan of the paper is as follows: in Section we will summarize the 



OPE approach, while in Section |2J] we describe a known mathematical tech- 
nique that appears to be useful to regularize and disentangle I.R. divergences 
of integrals, the Mellin transform. In Section ||] we apply the approach to the 



Ising model, (all the bulk computations are confined to Appendix A|) ; in Sec- 



tion 3.1 we use nonperturbative as well as numerical inputs to estimate the 



parameters left unfixed by the OPE method; found results are discussed in 



Section 3.2. Conclusions are given in Section 



2 All order I.R. finite formulas 



The goal of the method presented in |7J (see also [§, |9], [UJ and ]TT|, [L2|, [13[ 



for many preexisting ideas) is to obtain informations on the short distance 
behavior of a conformal field theory perturbed by relevant operators, i.e. 
when to the conformal action Soft is added a perturbation of the form 



AS = -JdxY,X B O lB (x), (2.1) 



where \ B (Oi B ) are bare couplings (operators) and < Xi < 2 (x, = xo { , 
xo being the scale dimension of operator O). This is achieved by expanding 
in powers of the corresponding renormalized couplings A* (Taylor expansion) 
the so-called Wilson coefficients C^ b (r, A) that enter in the Operator Product 
Expansion 

< $ a (r)$ b (0) > A ~ J2 C c ab (r, A) < $ c (0) > A (2.2) 
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for the complete theory ($ a are deformations of the conformal theory oper- 
ators, the suffix A refers to the complete theory correlators; the dependence 
on A* of the Wilson coefficients will be omitted in the following). 

Assuming the regularity of the Wilson coefficients in terms of the renor- 
malized couplings A* (minimality of the renormalization scheme), the validity 
of an action principle for the derivative with respect to A* and the asymp- 
totic convergence of OPE, all reasonably satisfied by the complete theory, an 
I.R. finite representation has been given for the n th derivative of C c ah with re- 
spect to the couplings evaluated at A* = 0, involving integrals of (eventually 
renormalized) conformal correlators. 

More in detail, by repeated differentiation with respect to couplings of 
Eq. fl2.2|) and by use of the (renormalized) action principle 

d X i <[■■■} >x= J d 2 x <{O i -- ■] > A (2.3) 

(square brackets [• ■ •] meaning renormalization) one can first obtain some in- 
tegral representation for the multiple derivatives of the Wilson coefficients in 
terms of integrated correlators of the complete theory. The limit A* — ► can 
then be taken safely (after regularizing the integrals by some generic I.R. 
cutoff function 0r(x) such lim^oo Qr(x) = 1) because OPE asymptotic 
convergence guarantees that the sum of all (singularly I.R. divergent) contri- 
butions is actually finite. In the limit in which all I.R. cut off are removed 
(R n , ■ ■ ■ R\ > R — > oo below) one also can show by dimensional arguments 
that the OPE series truncates to a simple sum, J2t below. As a result it was 
obtained in that 

jld ll ---d ln C b aia2 <^ b X 00 }> 
b 

= lim { Idxx ■ • • ldx n Q Rl {xi) ■ ••QR n (x n ) x 

R— too J J 

<[:O in :---:O h : ($ ai $ a2 - £ Cj^)**,] > 

b 

b 

< [: O in ■.■■■:O i2 : ^X^} > +p. + • ■ ■ 

- IA • -^iC fdx n e Rn (x n ) < [: O ln : ^XJ > +p.} (2.4) 
b J 

where di = p. means all nontrivial permutations over labels of indices i, 
Xri is an arbitrary operator localized outside the sphere of radius R' » R — > 
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oo, chosen to give adequate boundary conditions (and possibly including 
powers of R'), and J2ti as explained before, is restricted to operators $j, such 
Xb<(D- x ik ) h (D - x in ) - x Xr - 

While each term in the sum in the right hand side is I.R. divergent, (and 
consequently requires a regularization) the complete right hand side is I.R. fi- 
nite and independent on the choice of the cutoff function (what we presented 
above is an "I.R. renormalization" more than a simple regularization). The 
structure of the expression is simple: the first term on the right hand side 
is a naive (generalized) perturbative one, while the others, containing lower 
derivatives of Wilson coefficients, play the role of nonlocal I.R. counterterms, 
naturally induced by the theory itself. (A conjecture for the existence of this 
mechanism can be found in |0| , in the general context of quantum field the- 



ories. Also a rigorous proof of the conjecture within the MS renormalization 
scheme for perturbative superrenormalizable quantum field theories can be 
found in [|TJ.) 

If the cutoff function is rotationally invariant (Qr(t) = 0_r(|?"|)) only 
scalar operators contribute in all expressions above. Also, if by dimensional 
analysis one knows a priori that no powers R° can be obtained from the I.R. 
counterterms available in that model, it suffices then to compute the first 
(naive) term of right hand side in some chosen I.R. regularization, and to 
keep only the regular term of its asymptotic expansion when Ri — ► oo, the 
singularity of the naive term being automatically killed by the I.R. coun- 
tertermsQ (which cannot give in this case additional finite contributions). 
Roughly speaking, one expects that this picture is realized when the I.R. 
counterterms do not need U.V. renormalization (i.e. are locally integrable), 
because in this case the renormalization scale /i is absent and there is no other 
possibility to obtain an adimensional quantity from the log R terms that are 
expected in general from the R° I.R. counterterms (apart from miraculous 
cancellations of log R divergences). This shortcut can be applied in particular 
to the case we will consider here. 

We conclude by observing that the complete correlators are then obtained 
by combining the derived expressions for the Wilson coefficients with those for 
the required operator Vacuum Expectation Values (eventually parameterized 
by unknown but fixed and universal constants, see Section |3J] and discussion 
in0). 



lr This idea was already present in the first order computations of 
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2.1 Mellin transform and I.R. divergences 



Having in mind that all we need are asymptotic expansions of (multiple) 

integrals in the limit of m, = 4 > above, it is worth introducing here a 

natural method to deal with these expansions: the Mellin transform [] (see 



[31 1 for a nice introduction). 

Given a function I(m) locally integrable on (0, oo), of order m a when 
m — > and when m — > oo, one can introduce in the complex strip 
—a < s < (3 its Mellin transform 

Us) = / —m s I{m). (2.5) 
jo m 

If I(s) has a meromorphic extension in the complex plane on the left of 
Ims = —a one can derive from the inverse Mellin transform the asymptotic 
expansion, 

I{m) ~ Y, Res { m ~ S h s )) _ _ ( 2 - 6 ) 

i 

with ai = a < (%2 < ■ ■ ■ are the powers of m in the asymptotic expan- 
sion of I(m) (notice that multiple poles give logarithms of m). In a word: 
singularities of Mellin transform drive the asymptotic expansion. 

In our formulas we have always (possibly nested) integrals of the form 

I(m) = J d 2 z&(m\z\)g(z) (2.7) 

(we refer from here on to the rotation invariant cutoff; notice also the slight 
and obvious change of notations). It is easy to prove the following convolution 
theorem: 

I(s) = 0(s)G(l - s) (2.8) 

where 

G(l-s) = [ d 2 z\z\- s g(z) (2.9) 



is essentially the Mellin transform of g (up to angular integrals). Thus in 
the cases in which G is known the integrals are easily done. Notice that 
if the Mellin transform of g strictly speaking does not exist but the Mellin 
transform of I(m) is well defined, one can try to bypass the problem extending 
the integral I(s) to a more general one, 7(s; 5), (in which g is now extended 
to some adequate g$) such that it reduces to the original for some value 
5 = do, and such that convolution theorem can be applied: if 7(s; 5) can 



2 This technique was essential for the analog analysis in perturbative superrenormaliz- 
able quantum field theories Jl3) (and references therein). 
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be analytically continued to S = So, the Mellin transform of I is obtained 
accordingly. 

In particular if the I.R. counterterms give no finite contributions, (see 
Section |2[) then the derivatives of Wilson coefficients can be unambiguously 
obtained simply by taking the residue of the first, "naive perturbative" , term 
in 

In [Appendix A| we give the analytical expression of I(s) for a general case 



that allows us to compute all the integrals involved in the first order deriva- 
tives of Wilson coefficients of primary fields of the conformal field theory 
underlying the critical Ising model. 



3 Wilson Coefficients for the perturbed Ising Model 



The scaling limit of the lattice (D=2) Ising model (see |2(J and references 
therein) at the critical point, T = T c , H = is described by a continuous 
unitary conformal field theory, technically .M(3/4), The operator space 
of this conformal theory is generated by the primary operators l,a, £ of 
dimension x = 0, 1/8, 1. The related fusion rules are: 

[a][a} = [l} + [£], [£][£] = [1], [a}[£] = [a]; (3.1) 

implying 

< [£} k [l] m [a] 1 >= (3.2) 
for I odd. Also Kramers Wannier duality implies 

< [£]»[!]' >= (_)" < [£} n [l] 1 > . (3.3) 

The object of our study will be the extended theory in presence of a 
magnetic field: 

S = Si s - mg + J had 2 z. (3.4) 

Notice that with this notation each derivative dh will be realized by the 
insertion of the operator — / d 2 zo(z) (action principle). 

By application of the general formula ( |2.4| ) and of the selection rules ( |3.2| ) 
and (|3.3| ) we obtain the following (nontrivial) first order order relations 

- d h CZ a {r) < a a R , >= f d 2 z < a z a R ,[a r a Q - C^(r) - C^(r)£ ] > (3.5) 



- d h C% e (r) < ct (Tri >= J d 2 z < a z a R ,{£ r £ - C l ££ (r)} > 



(3.6) 
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- d h Cl £ {r) = fd 2 z < a z [a r £ - C^(r)<7 ] > (3.7) 

-d h C £ a£ {r) < £ £ R , >= J'd 2 z < a z £ a [a r £ Q -Cl £ {r)a Q -Cf £ {r)al] > (3.8) 

(where the prime on the integrals denotes a rotationally invariant I.R. regu- 
larization, all limits have been omitted and we defined cr 1 = L_ 1 L_ 1 a). 

It is easily realized by dimensional analysis that in all expressions above 
no finite contribution can arise from I.R. counterterms: we are thus in the 



case described in Section [2.1| , where computations are easier. 

To reconstruct the short distance corrections of complete correlators we 
must have the expressions of the VEV of the lowest dimensional operators, a 
and £. As well known in perturbative quantum field theories and from direct 



Renormalization Group considerations, see e.g. Hl9| , composite operators 
VEV's evolve by RG as 

| < [<& a ] >= It < [d> 6 ] > (3.9) 

in which I is the logarithm of the scale and the matrix T b a contains only pos- 
itive integer powers of the renormalized couplings A* and forbids the mixing 
with higher dimensional operators (it is lower diagonal if operators are or- 
dered with increasing dimensions). From these properties and from the fact 
that h has dimension 15/8 and that it is not renormalized it follows by simple 
dimensional considerations that a and £ do not mix with other operators, 
that T^, T £ reduce to the fixed point scale dimensions and that no logarithms 
appear in their VEV's. We can thus parameterize the VEV's as 

< a > h = 4r|/*| 1/1B Sign(/0 < £ >h= A £ \hf l \ (3.10) 

(where it has been also kept into account the fact that < a > h (< £ > h ) have 
to be odd (even) in h as easily seen directly from the corresponding lattice 
expressions) . 

Defining C = C | r | dime ^ p AB _ < j^q >ft | r |dimAB an j ^ sca ]i n g var j_ 
able t = \h\\r\ 15 / 8 we have the following expression for the short distance 
(small t) behavior: 

F aa = C^ a + A £ C§^ 15 + A a d^ a t w / 15 + 0(t 2 ) (3.11) 
F ee = C^ £ + A a d^ £ t 16 / 15 + 0(t 2 ) (3.12) 
KF a£ = A a Cf £ t l ^ + d^l £ t + A £ d^ £ t 2 ^ + 0{t^ (3.13) 
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where k = Sign(/i) and 



tic* - 1 ^(-i/2)r(7/ 4) \ 2 

^ Cctct " 32^ T(5/4) ) {6 - 1A) 

Wh = (3.15) 

^ 1 sin(37r/8) r 2 (ll/8) 

dfcC " £ " 4sin(157r/8)r 2 (15/8) ^ lb) 

fCr* ■ ( /srt ■ /Q J fr(-ll/8)r(7/8) V 

d h C£ £ = sm(7r/8) sm(37r/8) I — fTZw^ 1 ( 3 - 17 ) 



have been computed in |Appendix A| (see also [Appendix lj| for zeroth order 



expressions) and the only unknown constants, up to now are As and A a . 
While in general the OPE approach leaves some universal unfixed constants, 
in this case the (approximate) value of A a (As) can be obtained from the 
existing knowledge, see next Section. 



3.1 The missing VEV's 



The coefficient A a in (|3.10|) can be easily extracted by use of the results 
of |2l| (see also references therein), where the integrability of the A^(3/4) 
+ $12 (= er) was used essentially (in particular the Thermodynamic Bethe 
Ansatz and the method of external fields have been used). 



We report here Eqs. (2.9), (3.15) of ET|, 



M 



ATI 2 



4/15 



7 2 (3/16) 7 (l/4) 
4.404908579981566037 • • • h s/15 



4sin(7r/5)r(l/5) 8/15 

r(2/3)r(8/5) 



(3.18) 



M 



sin(7r/5) 



\2 sin(7r/5) J 8 sin(7r/3) sin(87r/15) 
-0.061728589822368- ■ - M 2 , 



(3.19) 



where M is the minimum mass of the particle of the theory, j(x) = T(x) /T(l— 
x), h > in this Section and the bulk energy e is defined by the partition 
function of the theory I.R. regulated in a cylinder of dimensions L, R: 



Z(h) =< e 



-Jha 



>~ e 



-eRL 



L,R 



oo. 



(3.20) 
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Noticing then that 



< a > h =< ae-J ha > / <e-f hu >~ Yim(-d h log Z(h)/LR) (3.21) 
it follows 

< a > h = d h e = A a h 1/15 (3.22) 

where 

A = 8 ( *** V /15 Sin ^ ( r(1/5) V(3 23) 
15 l v7 2 (3/16) 7 (l/4)y' sinf sinff Vr(2/3)r(8/15)y i ' ) 

= -1.27758227605119295- •• (3.24) 



On the other hand no nonperturbative informations are available (in our 
knowledge) for the VEV of operators different from the perturbation, o. 

In principle, numerical estimates of the constant As can be obtained by 
comparison with Monte Carlo simulations on the lattice for the normalized 
(connected) correlator < a(R)a(0) > c / < a > 2 in the regime in which 
R » 1 and t 8//15 oc R/£ is not too big, i.e., being £ oc H~ 8 / 15 , in a small 
field (R, £, H being lattice spacing, correlation length and magnetic field and 
the lattice spacing being assumed everywhere to be one). In practice with the 



available data, ||22|| , one has only a small window of applicability (£ ~ 15.5 for 
H = .001). To improve convergence and enlarge this window of applicability 
we fitted, together with the powers given in Eq.( |3.1ip also the next to leading 



terms t 2 ,t 32 ^ 15 (notice that logarithmic terms are absent for the same con- 
siderations given in Section |^). From the known first two coefficients in the 
small t expansion for the continuous expression of < a(R)a(0) > c / < a > 2 
we can fix the overall normalization (i.e. the lattice magnetization < a >i a t) 
and the correlation length £, related to continuous variables by R/£ = Mr 
(M being the mass in Eq. ( |3.18|) ). Notice that the predictions (|3.11|) - (|3.13|) 
refer to the complete correlators, and that the normalization factor between 
the continuous spin VEV and the lattice one (physical magnetization), is 
negative in our conventions, see also [g^l • 

At this point all is fixed and we get a prediction for the coefficient of 
£8/15 (Ag) and for the coefficient of t 16 / 15 (that actually is known exactly). In 
practice we selected only those fit that give a good prediction for the known 
coefficient, fixing in this way the above mentioned (but unknown) window of 
applicability: a good representative of which is (for H=.001) the R = 5 — 54 
range that corresponds to 

A s ~ 0.321 (3.25) 
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(see Eq. QB.lD for our normalization) with an error of order few percent (es- 
timated roughly from the variation of this quantity among different "good" 
fits). As a byproduct we also get 

f ~ 15.4 < a > lat ~ .634 (3.26) 

in agreement with the estimate £ = .38(l)if ~ 8 / 15 of [2^] and < a >\ a 



lat" 



1.003(2) of |H (with larger lattice). 

As an independent check that our estimate is correct we can use the 
known exact sum rule [IB], [T7|, ^4] 

A CT = - 2 ^ (2 ~ 2ACT) / <?r < a(r)a(0) > Kc . (3.27) 
47T < a > J 

We split the integral into two pieces, Mr < A and Mr > A. We estimated 
the first term by use of our expression ( |3.11 ) and the second one by use of 



the first three terms of the long distance expansion 

3 1 

< a(r)a(0) > ft , c ~ -{^fMamr) (3.28) 



where coefficients Cj and form factors F? can be found in [16[] (and references 
therein). By use of ( |3.27| ) we can give an expression of As in terms of A. 
Clearly for A small we have a big error from the approximated long distance 
expansion, while for small A we have a big error from the approximated 
small distance expansion. If (by a minimal sensitivity criterion) we choose 
the value of A that minimizes A? (A ~ 1.5) we obtain the value A? ~ .322, 
with an error of order 2% (estimated by repeating the procedure with only 



one mass term) that is compatible with Eq.(3.25). 



3.2 Comments 

A first glance at our results shows that a coefficient of the expansion for 
< ££ > ( |3.15|) vanishes nontrivially (unexpectedly from the point of view of 



the selection rules of the conformal field theory). This might sound strange, 
but one must bear in mind that the complete theory is an integrable model 
with powerful symmetries: we can thus interpret the found zero as a (short 
distance) signal of these symmetries and possibly of the existence of nontrivial 
differential equations satisfied by the correlator (analog of Painleve equations 
for spin-spin correlators in thermal perturbations of the Ising model, ) . 



Another important feature of our prediction is the presence of fractional 
powers of the coupling (or equivalently of the scaling variable t = \h\ |r| 8 / 15 ): 
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in the OPE approach these powers come out naturally from the VEV of 
operators (times some integer power of t coming from the Wilson coefficients, 
analytic in t). In particular the fractional correction t 8//15 that is found in 
( 3.1 H) for the < a a > correlator is expected if one introduces as infrared 



cutoff in perturbative expansion for the correlators the correlation length 
R c ~ |/i|~ 8 / 15 (see [|14], M). The existence of this fractional contribution to 
the spin spin correlator can be recovered from the results of |12| (based too 
on OPE and Callan-Symanzik equations in the 4 — e framework) by inserting 
the exact values of critical indices given by the conformal theory. 

Also it is interesting to see that the scaling of the < aS > correlator 
(vanishing at the critical point), is dominated by a fractional power, t 1 / 15 , 
coming from the VEV < a >. 

We must remark here that our results are in disagreement with an existing 
expression of the short distance behavior of spin spin correlator in the same 
context, PS| coming from the effort to apply in this case the perturbative 



expansion for correlators developed in || and to obtain the expected frac- 
tional powers of h as well. Main difference in with respect to our ( 3.1 1| ) 



is the presence of powers of logi in the first terms that, in our mind, are 
not motivated due to the absence of renormalization of lowest dimensional 
operators a, S (and of L_ 2 -^-2i as well, that is the next scalar operator that 
can have a nonzero VEV in this theory, as shown in [§, |i~9|), see discussion 
near Eq.( |3.9| ). Also the power t 24 / 15 present there is unmotivated from the 
OPE point of view (there is no corresponding VEV). 

In Figure 1 we give a comparison of the numerical data for < oo > with 
our prediction Eq.( 3.1l| ) (without the higher order terms added in previous 



Section to improve the fit) where the estimated values ( |3.25p as well as ( |3.26| ) 
are used. What appears immediately is that the convergence of the expansion 
is quite good: with the first order corrections to Wilson coefficients, we can 
achieve a relative error of order 2% at R/£ = Mr ~ 1. 

No numerical results exists in our knowledge for < aS > and < SS >. 
Nevertheless in Figure 2 we show a comparison between our (short distance) 
prediction and the (long distance) form factor result (up to the first eight 
terms), |T7|], in the case of < oS >. It appears from Figure 2 that there is a 
good agreement of two methods in the intermediate region 1 < Mr < 1.5 and 
a reasonable evidence of convergence of the long distance expansion towards 
our result in the region Mr < 1, as it is expected. This agreement could be 
regarded as an independent test of the form factor clustering property that 
have been assumed in |TJJ to get the form factors for < aS > and < SS >. 
(The origin of this property will be explained in ||24|| ). 
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4 Conclusions 



We computed the short distance behavior of all the correlators of primary 
fields of the critical two dimensional Ising model perturbed by a magnetic 
field up to 0(t 2 ) excluded, t = \h\\r\ 15 ^ 8 being the scaling variable. The 
I.R. finite OPE approach, extensively developed in |7j has been used. This 
method does not give informations about VEV's by itself. An unfixed con- 
stant for < a > has been obtained by use of nonperturbative results coming 
from the integrability of the model (TBA), while by comparison with the 
existing numerical data we got the nontrivial estimate (|3.25| ) of < £ > (not 
predicted from TBA). 



We point out once more that with the OPE approach the expected |14j 
fractional powers of the coupling h are naturally obtained from the operators 
VEV's, while in alternative perturbative direct estimates of the correlators 
(I.R. regularized in some way) these non integer powers are absent. In par- 
ticular we found immediately the expected t 8 ^ 15 correction to the spin spin 
correlator, |14|], and we predicted that the dominant term of < a£ > is the 



fractionary term t 1//15 (that should be observable in adequate lattice simula- 
tions). 

The result for < oo > when compared with the available data signal a 
good convergence of the approach: this encourages to search some resum- 
mation technique that could enlarge furthermore the convergence radius. 
Moreover the use of exact sum rules to extract informations from both short 



and long distance approximations that has been done in Section |3.1| , should 
stimulate the research in this direction, with the final goal of connecting the 
two regimes. 

We conclude by emphasizing that the OPE approach, being very general 
(integrability of the theory is not necessary to compute Wilson coefficients) 
and fastly convergent, can always be seen as a bridge towards the study of non 
integrable perturbations and as a test for the ansatz behind integrable theo- 
ries predictions (as done in Section |3.2j ). It is furthermore an open question 
if the powerful Coulomb gas technique of |2(| fZ7§ and the IR regularization 
of || could be married with the OPE approach to reach higher orders of the 
expansion in this and other statistical models. 
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Appendix A Computation of the integrals. 

First of all we give the general expression for the Mellin transform (see 
Eq.( [2.5D ) with respect to m (I.R. cutoff) of the integral: 

I a ,p, 7 (m; x) = J d 2 w&(m\w\)\w\ 2a \w - x^lw - 1\ 2/3 (A.l) 

in which m > 0, x G (0, 1), a, 7 > —1 (to have local integrability) and we 
will fix the I.R. cutoff function G(t) = e~* (such that its Mellin transform is 
0(s) = T(s)). By use of the convolution theorem (|2.8|) it easy to obtain 



4a7(«; x ) = F(s)D(a - s/2, p, 7, x) (A.2) 

where the integral 

D(a,b,c,x) = [ d 2 w\w\ 2a \w - x\ 2c \w - l\ 2b (A.3) 



has been computed in [E7| by use of contour deformation in the variable Imw 
and found to give: 

r^/ , \ S(a)S(c) . T l2 S(b)S(a + b + c) . T l2 , . 

D(a,b,c,x) = ) v f /ox 2 + y T \ - ^ioo 2 A.4 

o(a + c) D(a + c) 

where 

5(a) = sin 7ra 
Jo- = I* duu a (x-u) c (l -u) b 





x 1+ « +c ^±^±^F(-6, a + 1; a + c + 2; x) 
L{a + c + 2) 



Jioo = J duu a (u — x) c (u — 1) 

T(-a-b-c- l)T(b+ 1) 
T(— a — c) 

and J 1 = 2J1 is the Hypergeometric function (see e.g. |^[)- From the 
Mellin transform ( |A.2|) we can obtain the asymptotic expansion of the original 



F(—a — b — c — 1, — c; —a — c; x) 



integral (|A.1|) when m — > 0, by use of 



13 



All the required integrals in (|3.5| - |3.8|) can be then computed by use of 
the general expression ( A.2|) . See [Appendix B| for the conformal correlators. 
In our particular case we know that the I.R. counterterms cannot give m° 
corrections, see Section |3| so that the derivative of the Wilson coefficient is 
obtained by extracting the m° contribution of the naive perturbative term 
(the first of right hand side) in ( |3.5| - |3~8| ), i.e. by computing the residue of its 
Mellin transform in s = 0. Notice that due to the absence of m°logm I.R. 
counterterms, there cannot be log m factors in the naive term and correspon- 
dently its Mellin transform will have only simple poles. While (|3.7| - |3.8|) are 
trivially reconduced to ( |A.2| ), some few words must be spent for ( |3.5| ), i.e. for 
the Mellin transform of four spin correlator ( |B.6| ). In the limit R'm —>■ oo (R' 
being the argument of the boundary operator <jri) we have in the correlator 
the simplification x = ( p,Zy\ = - (r = 1) so that it is better to change 



variable and obtain (in notation of the convolution theorem (2.5)): 



G(l-s) 



d x\x\ 



|1 — x\ 



-26 



1 + Vl 



X 



x 



|). (A.5) 



As explained in Section [2.1| we have introduced an additional parameter 
5 to justify the use of convolution theorem. The analytic continuation to 
5 = 1/8 will give the wanted Mellin transform of the four spin correlator 
with exponential cutoff. The final step to reduce (|A.5| ) to (|A.2 ) is the simple 
substitution w 
for w 



t/T— x and the observation that the integrand is invariant 
—w. By use of these considerations Eqs.( 3~T4] - |3.17 ) follow. 

As a general check of the regularization independence our approach, we 
report the alternative derivation of dhC° a obtained regularizing the integrals 
by restricting them to \z\ < R. We will not use the Mellin transform tech- 
nique and we will keep explicitly the I.R. counterterms to show explicitly 
how cancellations work. 

as 



In the limit R! j R — ► oo (being x = 1/z) we can rewrite ( p75| 



- d h C%„ 



d 2 x , 



|x|>l/i? \X\ 



\9 W 



+ \9 



1, 



\x 



(A.6) 



where we omitted the remaining overall limit R — > oo in left hand side and 
defined 



9 ± (x) 



\ 2(1 - x y/*' 



(A.7) 



Unfortunately we have not a closed expression for the integral ( |A.6[ ). 
However the calculation can be performed by series, splitting the integral in 
two regions, 4 < \x\ < 1 and \x\ > 1. The contribution of the first region is 
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then obtained expanding 



oo 



//•i.r) X /// g-(x) = £<?-^ +1/2 (A.8) 
i=o j=o 

(note also that (g±(x))* = g±(x*)), exchanging the series with the integral 
and performing simple integrals of powers of x x*. At this stage the role of 
I.R. counterterms is essential to give a finite R — > oo limit. The contribution 
of the second region is obtained after making the change of variable x = — , 
from which we obtain 

2 / d 2 w\wf 2 \h{w)\ 2 (A.9) 

J\w\<l 



h{w) 



\ (TT^r (A ' l0) 



for the sum of g ± contributions (integrals of I.R. counterterms are easily 
performed). Again developing 

oo 

h(w) = Y,hjW i (A. II) 

3=0 

the integral can be easily performed by series. 
The final result is 



- i I v-/ l^l 2 \9j \ 2 \ 3 4 l x . 2|/i,-| 



(A.12) 

where the coefficients /ij satisfies the recursive relations 



(j + l)(j + ±)g+ +1 - 2j(j - ±)g+ + (j 2 - lj + t = 

3 13 3 5 

(j + 1)0' + gtel - 2(j + -)(j + g)^7 + (f - - a j + ^)<77_! = 

(j + l)(j + 2)h j+2 +j(2j - l -)h 3 + (j 2 - + ^)V 2 = (A.13) 
due to the existence of the following differential equations satisfied by g ± , h: 
(x 3 - 2a; 2 + x)g ± + {-x 2 - -x + -)g ± - —xg ± = 
(1 + 2u> 2 + w A )h" + ^{w 3 + w)h' - ^-w 2 h = 0. (A.14) 
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As they stand the series are very slowly convergent, (g^, hj 



-7/8 



) SO 



to increase the convergence of the series we added and subtract to the j th 
term its asymptotic powerlike behavioif] that can be resummed and gives 
well known Lerch functions ($(z,s,a) = J2 n +a^o ( n +a) s ' n = ^> ^' ' ' ") 
are obtained with high precision from Mathematica. 

After numerical resummation we obtain then 

d^CZa = -403745631213448123 ■ • • (A. 15) 

that agrees up to the obtained precision with the analytic result (|3.14Q. 



Appendix B Conformal correlators 

To fix notations we give here the expressions for the fixed point (confor- 
mal theory) quantities involved in our computations. See ] 
derivation. 



iT| for their 



Wilson Coefficients: 



Cl(r) 



Ch{r) - 
1 



47r 2 |r| 



CL(r)=rr\r\^ 
1 



47r|r| 



Useful correlators: 



< a(z 1 )a(z 2 )£{z 3 ) >- 



1 


\Z12 


3/4 


4tt 


1*13 1 


^23 



(B.l) 
(B.2) 
(B.3) 

(B.4) 



< a{z 1 )a{z 2 )£ {z 3 )£(z 4 ) > 





^12(^32 + Z42) - 


- 2^32^42 2 


16tt 2 


^42^32^41^31 


^43 


2 


Zl2 


1/4 



;b.5) 



1 4" VT" 

< a(z 1 )a(z2)cr(z 3 )(T(z 4 ) >= | (1 - x)z 12 z u \~ 1/4 (\ — 



X , 



where x = z 12 z u / (z 13 z 24 ), % = Zi - Zj. 



(B.6) 



3 The first two terms are obtained from the contribution of the singularities nearest to 
zero (Darboux theorem, see [BOl M) while the others can be conveniently extracted by 



expanding in j 7 / 8 + fe / 2 the recursive relations (A. 13) up to the desired order k. 
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Figure Caption: 

Figure 1 Values of < aa > from lattice simulation (dots) and as predicted 
from OPE approach (solid line) versus R (lattice position) for H = .001 
(£-15.5). 

Figure £ Estimates of versus Mr: comparison between the OPE 

approach (solid line) and the form factor method with one (dashed line) three 
(dot-dashed line) and eight (dots) terms. 
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